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We present the results of a search for the radiative decay B° J/i'l in 8- data set containing 
123 million T(4S') BB decays, collected by the BABAR detector at the PEP-II asymmetric- 
energy e'^ e~ storage ring at SLAC. We find no evidence for a signal and place an upper limit of 
B{B° J/ip-)) < 1.6 X 10"'' at 90% confidence level. 

PACS numbers: 13.25.Hw 



Rare decays are sensitive probes of possible new 
physics effects beyond the Standard Model. The de- 
cay i?" J 1^1^ is a very rare decay, with a predicted 
branching fraction of 7.65 x 10~^ The dominant 

mechanism is the exchange of a VF-boson and the ra- 
diation of a photon from the light quark of the B me- 
son (Fig. Possible new physics enhancements of the 
Jhh'y decay rate include a right-handed charged 
current llj or non-spectator intrinsic charm in the B'^ 
meson No prior search has been conducted for this 
decay mode. 

The data used in this analysis were collected with the 
JBABj4R detector at the PEP-II asymmetric-energy e~^e~ 
storage ring. The data sample contains a time-integrated 
luminosity of 113 fb~^ of "on-peak" data taken at the 
T{4:S) resonance, corresponding to a center-of-mass en- 
ergy of = 10.58 GeV, as weh as 12 fb"^ of "off-peak" 
data recorded at about 40MeV below this energy. The 
energy asymmetry between the low-energy positron beam 
(3.1 GeV) and the high-energy electron beam (9.0 GeV) 
produces a Lorentz boost of /? = 0.56 of the center-of- 
mass frame with respect to the laboratory frame. In 
this report, measurements are presented in the labora- 
tory frame unless otherwise specified. The z-axis points 
along the direction of the high-energy electron beam; po- 
lar angles (9) are measured with respect to this axis. 

The BABAR detector is described elsewhere here 
we provide a brief overview. The detector consists of 
five subdetectors. The tracking system includes a 40- 
layer, helium-based drift chamber (DCH) as the main 




FIG. 1: Feynman diagram of the leading-order contribution 
to B° ^ J/ip-y. 



[1] Charge conjugation is implied throughout this paper. 



tracking chamber, and a five-layer silicon vertex tracker 
(SVT) for precise reconstruction of track angles and B- 
decay vertices. The tracking system covers the polar an- 
gular region 0.41 < 9 < 2.54 rad (86% of the solid angle 
in the center-of-mass frame). Charged-particle identifi- 
cation, particularly K/t: separation, is provided by the 
DIRC, a ring-imaging Cherenkov detector. Electrons and 
photons with energy greater than 30 MeV and a polar an- 
gle within 0.41 < 9 < 2.41 rad (84% of the soHd angle in 
the center-of-mass frame) are detected by an electromag- 
netic calorimeter (EMC) with energy resolution at 1 GeV 
of 2.6%. These four subdetectors are contained inside a 
magnetic solenoid, which supplies a 1.5-T magnetic field 
for tracking. The fifth subdetector, a segmented iron flux 
return (the IFR) , surrounds the magnetic solenoid and is 
instrumented with resistive plate chambers for muon and 
K1 identification. 

Particle candidates are either charged tracks in the 
tracking devices (SVT, DCH), or "clusters" — groups of 
adjacent hits — in the EMC or the IFR. Each charged 
track is tested to see if it comes from the same parti- 
cle as one of the clusters, and if so, it is matched with 
that cluster. Charged-particle candidates are thus either 
stand-alone charged tracks or track-cluster pairs, while 
neutral-particle candidates are clusters not matched with 
charged tracks. 

The analysis proceeds as follows. We use simulated 
signal and background samples to derive an optimized 
set of selection criteria, and to estimate the fractions of 
signal and background events that pass the criteria. We 
then apply the selection criteria to the data sample and 
calculate an upper limit on the branching fraction for 
B'^ — !• J/tpj. The simulated signal sample contains only 
B'^ J/i^^ events in which the J/ij} meson decays in 
the J/il) t^l~ mode, where I denotes an electron or 
a muon. The sources of background in this analysis in- 
clude background from B decays and from continuum 
quark production (e+e" qq, where q — u, d, s, c), the 
latter being three times the size of the former. The sim- 
ulated background sample contains both of these types 
of background. 

To obtain a i?i?-enriched sample, we impose require- 
ments optimized independently of this analysis and used 
in many other i3-decay studies. Events are required to 
have visible energy greater than 4.5 GeV and a ratio of 
the second to the zeroth Fox- Wolfram moment 0, i?2, 
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less than 0.5. We reconstruct a primary event vertex 
from the charged tracks and require that it be located 
within 6 cm of the beam spot in the direction parallel to 
the beam line, and within a transverse distance of 0.5 cm 
from the beam line. The beam spot RMS size is approx- 
imately 0.9 cm in z, 120 /im horizontally, and 5.6 /zm ver- 
tically. There must be at least three tracks in the fiducial 
volume satisfying the following criteria: they must have 
transverse momentum greater than O.lGeV/c, momen- 
tum smaller than lOGeV/c, and at least 12 hits in the 
DCH; and they must approach within 10 cm of the beam 
spot in z and within 1.5 cm of the beam line. Studies 
with simulated samples indicate that these criteria are 
satisfied by 96% of T{^S) BB decays. 

Candidate B^ J/i^l decays are reconstructed as 
follows. A candidate is formed from a J/ij} and a 
photon candidate. The J/il) candidate is reconstructed 
in the low-background, high-efficiency J/ij] t' mode 
only. Electron candidates are identified using the ratio of 
calorimeter energy to track momentum [E /p)^ the ioniza- 
tion loss in the tracking system (dS/dx), and the shape 
of the shower in the calorimeter. Whenever possible, pho- 
tons radiated by an electron traversing material prior to 
the DCH (0.04 radiation lengths at normal incidence) are 
combined with the track. These bremsstrahlung-photon 
candidates are characterized by an EMC energy greater 
than 30MeV and a polar angle within 35 mrad of the 
electron direction, as well as an azimuthal angle either 
within 50 mrad of the electron direction, or between the 
electron direction at the origin and the azimuth of the 
impact point in the EMC. Muons are identified by the 
energy deposited in the EMC, the compatibility of the 
track formed by the hits in the IFR with the extrapola- 
tion of a track measured in the DCH, and the amount 
of iron penetrated by this track. Studies of data-derived 
control samples show that at a typical lepton momentum 
of 2 GeV/c, the efficiency of the electron (muon) identifi- 
cation criteria is 93% (83%), with a pion misidentification 
probability of 0.2% (8%). Photons are neutral candidates 
with characteristic electromagnetic shower shapes in the 
EMC. To determine the photon direction, we assume that 
the photon candidate originates at the J/'ip £^£^ ver- 
tex. 

We use the simulated samples to derive an optimized 
set of selection criteria for 5° J/ipl events. For the 
optimization we minimize the ratio ^yeb/ss, where Sb and 
Es are the respective efficiencies for background and sig- 
nal events to pass the selection criteria. The optimized 
selection criteria are described below and summarized in 
Table D 

To identify and select B candidates we use the kine- 
matic variables AE and toes- The energy difference 
AE is given by AE = {2qr ■ Qb — •s)/2y^, where qr — 
{Er,pr) is the four-momentum of the T(4iS') as deter- 
mined from beam parameters, qs = qj/ip +97 = {Eb,Pb) 
is the reconstructed four-momentum of the B candi- 



date, and s = is the squared center-of-mass en- 
ergy. The energy-substituted mass toes is given by 
Toes = v/(s/2 +pr ■ Pb)^/E^ - bsP- The advantage of 
using AE and toes to impose the kinematic constraints 
for B decays is that these quantities are largely uncor- 
related and make maximum use of the well-determined 
beam four-momentum. For the optimization and back- 
ground studies we use only events that fall within the 
"analysis window" defined by 5.2 < toes < 5.3GeV/c^ 
and |Aii^| < 0.30 GeV; this defines the range of the his- 
tograms in Fig. 13 A perfectly reconstructed B^ J/ipl 
decay should have AE = and toes = "^b ■ Therefore we 
demand that 5° J/'4'l candidates fall within the "sig- 
nal region" in the AE vs. toes plane defined by 5.270 < 
Toes < 5.290 GeV/c2 and -0.05 < A£; < 0.08 GeV. In 
Fig. 12 the signal region is indicated by a box. 

We reject continuum background using a number of 
topological variables to distinguish between continuum 
events, which tend to be highly directional, and i3-decay 
events, which tend to be spherically symmetric. We de- 
termine the thrust and sphericity axes of the particles not 
used to reconstruct the B candidate, and demand that 
the angle 9t (Osph) between the thrust (sphericity) axis of 
these particles and the thrust (sphericity) axis of the B 
candidate satisfy |cos0f| < 0.75 {\cos9sph\ < 0.85). We 
further demand that the polar angle 9b, the angle be- 
tween the beam direction and the flight direction of the 
B candidate in the e~^e~ center-of-mass frame, satisfy 
I cos 6*3 1 < 0.90. This demand accepts or rejects events 
based on the angular momentum constraints of the de- 
cay T{AS) BB. Finally, we also tighten the R2 re- 
quirement to i?2 < 0.45. Studies both of simulated back- 
ground and off-peak data indicate that the fraction of 
continuum events satisfying these criteria is negligible. 

We reject background from B decays using J/iJj and 
photon selection criteria. For the J/tp selection, the 
invariant mass of the £^£~ pair of the reconstructed 
J/'ip £~^£~ decay is required to fall close to that of the 
known J/ip mass 3.06 < m{e+e-) < 3.12GeV/c2 
for J/i/j e+e~ candidates and 3.07 < m{fi^^~) < 
3.13GeV/c^ for J/i{} ^ fi~^^~ candidates. We require 
that photon candidates satisfy LAT < 0.35, where LAT 
't'I is a shower-shape variable used to distinguish be- 
tween electromagnetic and hadronic showers. In addi- 
tion, we constrain the photon direction to the region 
cos^^ > —0.35. 

The main source of photons in BABAR is the decay 
of neutral pions, so we apply a veto to reject photons 
from TT*' — 77 decays. We reject events in which the 
B^ — !■ J/V'7 photon candidate combined with any other 
photon candidate forms a pair with an invariant mass 
within 20 MeV/c^ of the neutral pion mass . The signal 
efficiency of the optimized selection is estimated from the 
simulations to be = 0.102 ± 0.010. 

Of interest in the background studies are the events 
that pass all of the selection criteria except for the re- 
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TABLE I: The selection criteria. 



Variable 



J/tp mass 



Requirement 

3.06 < m{e+e-) < 3.12GeV/c^ 

3.07 < m(/i+/i-) < 3.13GeV/c2 
LAT < 0.35 
cos^T, > —0.35 



photon LAT 
photon angle 

7r° veto reject 0.115 < m^pair 
Fox Wolfram moment R2 < 0.45 

thrust angle |cosSt| < 0.75 

sphericity angle \cos9sph \ < 0.85 

B polar angle | cos 9b \ < 0.90 
signal region 



< 0.155 GeV/c 



5.270 < mES < 5.290 GeV/c^ 
-0.05 < AE < 0.08 GeV 



quirement to fall within the signal region (Fig.|2t). Most 
of this background is concentrated in the low-Ai? region 
of the A_E-mEs plane. The asymmetry of the signal re- 
gion in AE ensures that the majority of these events fall 
outside of the signal region. The small fraction of this 
background in the signal region is due primarily to B'^ 
— > J/ip 7r° decays in which a photon from tt" — s- 77 is 
misidentified as a B'^ J/ip-j photon. This usually oc- 
curs when the other photon in the reconstruction falls 
below the 30MeV energy threshold. There is also back- 
ground from i?° J/^ decays, due to K° 37r° 
decays in the EMC for which the six resulting showers 
overlap and are incorrectly interpreted as a shower from 
a single photon. 

We estimate the background using a large simulated 
sample distinct from that used to optimize the selection 
criteria. Each event in this sample contains either a i? — )■ 
J/ipTT^ or a B ^ J/ipKl decay. After normalizing to the 
data luminosity we obtain background estimates of 0.59 
in the B J/tp-K^ mode and 0.12 in the S -> J/tpK^ 
mode, resulting in a total background estimate of rib = 
0.71 ± 0.31 events. The contributions to the uncertainty 
are discussed below. 

To validate the simulated-background modeling we 
perform several cross-checks. We compare background 
estimates from simulations and from on-peak data, out- 
side the signal region but in the analysis window. The 
results are consistent both when the estimates are ob- 
tained with all of the selection criteria applied, and when 
the estimates are obtained with all of the criteria applied 
except for the pion veto. In addition, we compare the 
background estimates from off-peak data and from simu- 
lated continuum background in the full analysis window. 
In both cases, no events pass the selection criteria. 

The relative systematic errors in the signal efficiency 
and in the background estimate are presented in Ta- 
ble ^ For both Eg and rib there is statistical uncer- 
tainty in the number of events passing the selection. 
The uncertainty in the background estimate also in- 
cludes uncertainty from the number of T(4S') in the 



TABLE II: Summary of relative systematic uncertainties on 
the signal efficiency Ss and background estimate nt. 

Uncertainty (%) 



Source 


£s 


rib 


AE correction 


7.6 


33 


Tracking 


4.0 


3.4 


Lepton ID 


3.5 


3.2 


Neutral ID 


2.5 


2.1 


Statistics (simulated samples) 


1.5 


24 


BiB° ^ J/V-ttO), B{B° ^ J/tpKl) 


N/A 


15 


B{J/iP £+£-), and B[r{4S) B°B°] 






Total 


9.8 


44 
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FIG. 2: A_E-mES distributions in the analysis window for (a) 
on-peak data, (b) simulated signal, and (c) simulated back- 
ground. The signal region is indicated by a box. The sample 
in (c) is the original background sample used to optimize the 
selection; the equivalent luminosity of this sample is about 
nine times that of the 113 fb~^ sample in (a). 



certainty in the following branching fractions. B{B^ 
J/ipTT°) and B{B" J/il^Kl) are obtained from Ref. 0. 
B{J/ip ^e+£-) = 0.1181 ± 0.0020 is the sum of the 
J/V' e+e" and J/?/' ^ M^/^^ branching fractions 0| 
assuming fully correlated uncertainties. B[T{AS) 
B°B°] = 0.499 ± 0.012 is determined from Ref. as- 
suming that the 7(45*) decays 100% to BB. 

In addition, we correct for differences between sim- 
ulations and data, and each of these corrections con- 
tributes to the systematic uncertainty. The required cor- 
rections for tracking, lepton-identification, and photon- 
reconstruction efficiencies are derived from independent 
studies comparing the results from simulations with those 
from data control samples. Also, comparison of the AE 
distribution of K*'^j decays in real and simulated 

samples reveals a difference of about 28 MeV in the cen- 
tral value for AE between data and Monte Carlo. This 
effect is due to imperfect simulation of photon energy loss 
in the detector. B'^ J/ipl is topologically similar to 
B'^ —^ K*^j but has a lower photon energy, so we apply 
a correction of (22 ± 10 MeV) to AE in the simulated 
samples. As shown in Table HTl this AE correction leads 
to the largest systematic error in both the efficiency and 
the background calculation. 

No events in the signal region satisfy the final selection 
criteria (Fig.[2Ji). The probability of observing events 
when expecting a background of 0.71 events is 49%. In 
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the analysis window we observe 10 events in data, con- 
sistent at the 8% level with the expected background of 
5.7+/-1.0 events. 

We determine the upper limit on the branching fraction 
B{B^ —> J/4'l) by performing a Bayesian analysis with 
a uniform prior above zero. We define the likelihood for 
B{B^ J/ipj) as the probability that exactly zero events 
pass the selection, given that the mean expected number 
of observed events is 

^, = 716 + Nbo Ss B{J/^ ^ e+£-) ^(5° ^ J/^A7), (1) 

where Ngo = 2 Ny^is) S[r(45) ^ B°B°] is the number 
of B^ mesons in the data set. The analysis takes into 
account the imccrtaintics in Sg and rifc. The 90% con- 
fidence level upper limit, defined as the branching frac- 
tion value that separates the lower 90% of the area un- 
der the likelihood function curve from the upper 10%, is 
B{B° J/V't) < 1.6 X 10-'3. This limit is dominated by 
statistical errors; in the absence of systematic errors, it 
would improve by less than 0.1 x 10~^. 
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